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Abstract 

We study the loop induced production of charged Higgs boson associated with 
W-boson via photon-photon collisions at linear colliders (LO). The investigation is 
carried out under the two-Higgs-doublet model(THDM) and the minimal supersym- 
metric model (MSSM) using full on-mass-shell renormalization scheme. The numerical 
analysis of their production rates is presented with some typical parameter sets. The 
results show that the total cross sections in frame of the MSSM at linear colliders(LO) 
can reach 0.1 femto barn quantitatively. The contributions from supersymmetric sector 
can lead to an order of magnitude enhancement of the cross section in some parame- 
ter space. We conclude that the associated W^H^ production rate in photon-photon 
collision mode is comparable to that in e^e'' collision mode. 
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1 Introduction 



As we know the Higgs sector is one of the essential aspects of the Standard Model(SM) 
responsible for spontaneous symmetry breaking and the generation of masses for the 
fermions and gauge bosons. But until now the origin of electroweak symmetry breaking 
has not yet been directly tested by experiments. The extensions of the SM can lead to ad- 
ditional physical Higgs bosons. Particularly the general two-Higgs-doublet model(THDM) 
and the minimal supersymmetric model(MSSM) necessarily involve charged Higgs 

bosons H^. Because of its electrical charge, the charged Higgs boson is noticeably different 
from neutral Higgs particles. Therefore, experimental discovery of Higgs bosons including 
charged Higgs bosons is one of the important goals of the present and future colliders. 

Here we review briefly the production mechanisms for charged Higgs bosons at hadron 
colliders for a comparison between the production mechanisms at LC and at the TeV 
energy scale hadron colliders, At hadron colliders, there are several mechanisms which can 
produce charged Higgs bosons p. 

(1) The charged Higgs boson pair can be produced by qq H^H^ [|6| and gluon-gluon 
fusions at one-loop gg H^H^ at the LHC. Since the heavy ff^-bosons decay 
dominantly into quark pairs, the pair production process is always accompanied by serious 
QCD backgrounds. 

(2) When mjj± < mt — mt,, the bosons can abundantly be produced in decays of 
top (anti-top) quarks ti ^ bH^(bH^) from the parent production channel pp — > tt. The 
dominant decay channels in this mass range are H^{H^) f z^T-(r^'T-). 

(3) When > mt — mb, the single charged Higgs boson production can be via gb{gb) 
tH-{tH+) @, 55 ^ tbH-{tbH+)§, and qb{qb) q'bH+{q'bH+) processes Q. The 
sequential decay tb is known as a preferred channel for boson search. But 
signals in these processes appear together with large QCD background. The associated 
heavy production with W boson is another important channel in production of the 



heavy bosons which has been investigated in Refs. |11]|12]|13]. In this case the 
M^^-boson's leptonic decay may be used as a spectacular trigger. 

In the future linear colliders (LC) there are several mechanisms which can produce 
charged Higgs bosons: 
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[1) Pair production via e~^e H^H and e^e — > 77 — > H^H . The former has 



aheady been studied at the tree-level in Ref.|14] and one-loop order in Refs. |15]|16]. 



The latter production process is through laser back-scattered 77 collisions and has been 



studied in Ref.|17|. It is found that these two processes are prominent in discovery of the 
charged Higgs bosons. 

(2) Single production of charged Higgs boson in association with W gauge boson at 
e^e~ colliders provides an attractive alternative in searching for , which is kinemati- 



cally favored when 'mfj± exceeds niw- In Ref.[18| the associated production mechanism 
via direct e^e~ collisions has been studied in the THDM. Since there is no tree-level 
H^W^V couplings {V = ^ and Z^), both the H^W^ assiociated production processes 
via 77 and e^e~ collisions have no tree-level contribution in the THDM and the MSSM 
and occur at one- loop level in the lowest order. 

In this paper we concentrated ourselves on the investigation of the process e"'"e~ 
77 W^H^ including the complete contributions in a non-SUSY THDM and the MSSM 
at LC. Since both the W^H^ associated production processes via e^e~ and laser back- 
scattered 77 collision modes are all one-loop induced at the lowest order, and from the 
technical point of view there is no problem in realizing the luminosity of laser back- 
scattered photon beam approching to that of electron beam, the production rate through 
77 collisions could be competitive with the e^e~ colliding mode at LC. We arrange this 
paper as follows. In Sec. 2 we present the analytical calculation. In Sec. 3 we give some 
numerical presentations and discuss the numerical results of the processes e^e~ ^ 77 ^ 
W^H"^ . In frame of the MSSM we do numerical calculation in the minimal supergravity 
(mSUGRA) scenario The conclusions are contained in Sec. IV. Finally some notations 
used in this paper and the explicit expressions of the self-energies are collected in the 
Appendix. 

2 The Calculation of e+e" ^ 77 ^ W+H 

In this section we investigate the process 

e+e" ^ 77 ^ W^H~ 
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in frames of a CP conserving non-SUSY THDM and MSSM. In these models the cross 
section of e+e~ 77 W^H^ coincides with the process e+e~ ^ 77 ^ W~^H~ because 
of the charge conjugation invariance. Our calculation in this section is concentrated on the 
process e+e^ 77 ^ W^H^ , unless otherwise stated. But the numerical results of the 
total cross section in Section III involve both two processes i.e. e+e~ ^ 77 ^ W^H^. 
Hence our numerical results of total cross section contain a factor of 2 in contrast to only 
W~^H~ production. 

In our calculation, we adopt the 't Hooft-Feynman gauge, and use dimension regu- 
larization scheme for the calculations in frame of the non-SUSY THDM and dimension 
reduction scheme in frame of the MSSM. We denote the reaction of W^H~ production 
via photon-photon collisions as: 

where pi , p2 and k\ , k2 represent the four momenta of the incoming photons and outgoing 
and , respectively. The Mandelstam variables are defined as s = {pi +^2)^, 
i = {Pi — ki)^ and u = {pi — ^2)^- 

We depict the generic Feynman diagrams contributing to the subprocess 77 — W~H'^ 
in Fig.l in the frame of the MSSM. As shown in the figure, this is a one-loop induced 
subprocess. The diagrams with internal supersymmetric particles q, %+ and x° do not 
appear in the non-SUSY THDM. In Fig.l there are self-energy, triangle, box and quartic 
one-loop diagrams as well as the counter-term for the 7 — W'^ — H" vertex appearing 
in t-channel. The corresponding u-channel diagrams are not shown in this figure. The 
existance of vertex counterterm is due to the fact that this vertex is absent at the tree 
level, but is UV-divergent at loop contributions. The 7 — — H~ vertex counter-term 
can be visualized in Fig.l(b.l), (b.2). As this subprocess involves non-diagonal gauge 
propagators, the loop diagrams in Fig.l(a.l-a.l9) and the one-loop diagrams in Fig.l(b.3- 
b.5) are not UV-finite by themselves. Both the mixing H'^ — W'^ and H^ — G^ propagators 
have to be renormalized to obtain finite physical results. The explicit self-energy one-loop 
diagrams for mixings of — W'^ and — are plotted in Fig.2. Since the H'^ — 
self-energy diagrams on the external W boson have no contribution as a consequence of 
on-mass-shell W gauge boson, we do not plot them on Fig.l. By using the renormalization 
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condition as described in Refs. ||2^ the sum of Higgs tadpole diagrams and tadpole 
counter-terms vanishes. Therefore, there is no need to consider the tadpole diagrams in 
the self-energy calculations. 

The Higgs sector of a non-SUSY THDM or the MSSM consists of two scalar doublets. 
In reference [^] the complete set of Feynman rules for the MSSM is given. Here we follow 
their notations and conventions. We denote the two Higgs doublets as 



^2 



Theoretically, the self-interactions among charged Higgs bosons of two-Higgs-doublet 
extension models and the gauge bosons are completely determined by local gauge invari- 
ance. The kinetic terms read: 



■'kin 



Yu V ^ Yu 



{d^ + igi^B^ + ig2T''WpH2. (2.1) 



2 "'^ ' 

where are the SU (2) generators. and denote the U{1)y and SU (2) l gauge fields 
respectively, and gi and g2 correspond to their relevant coupling constants. Y^-. {i = 1,2) 
are the hypercharges of Higgs fields. Our one-loop calculation involves — and 
— mixing self-energies and their relevant counterterms. We adopt the on-mass- 
shell renormalization scheme suggested in Ref . [^] for the extended Higgs doublets. We 
define that the Higgs fields and vacuum expectation values are renormalized as below: 

Hi zH^Hi 

Vi ^ Z^j^^iv, + bvi) = (^1 + + -bZu^^ Vi, (2.2) 



hence 



Stanf) 6v2 dvi 1 



tan/3 V2 vi 2 
We replace of the free parameters and the fields in the kinetic terms of the Lagrangian 
Eq.(2.1) by the renormalized parameters and fields to obtain corresponding counterterms 
at the one-loop order in the Lagrangian, then we can read out all the relevant counterterms 
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for our process. We define the counterterms SZ^g ^^HW^ wliich behave as if they 
were the renormalized constants introduced as: 



(2.3) 



and 



Zf)V2w±±,^5,i/±, (2.4) 



where Z^ = 1 + SZ2^ is the usual SU{2)l gauge triplet renormalization constant. Then 
we get the relations 



-{5Zh2 - SZhi) =F ^ 



2' '''' ' tan/3 



5Zhg = 6Z^^l + 5Z^^l = ^-^{6ZH,-dZH,), 

6Zhw = Sz'j^l. (2.5) 

And the counterterm for the jWH vertex appearing in Fig.l(b.l,b.2) has the form as 

SV^WH = ieg^,ymwSZ\jlj, (2.6) 

which is obtained with the above replacement of the Lagragian terms C = e mwA^W~G^+ 
h.c. + ... in Eg. (PI. 

Denoting the bare self-energy as S, the renormalized one is defined as S = S — JS. 
Here we define the renormalized two-point Green functions for the — W'^ and — 
mixings as 

pz _ mjjj^ p-^ — 

G^^ = 5 — C^HG — SY^Hg)^ 9-, 

— — 

respectively. The counterterm 6T, can be derived from the bare self-energy expression by 

using on-mass-shell scheme. With an on-mass-shell charged Higgs boson, the real part of 

the renormalized mixing self-energy T,Hw\p'2=mjj± should be zero, i.e. 

ReT^Hw\rfi=n? = =^ oZhw = 9 • v^-''') 



For the renormalized self-energies Tjhw and TihGi we have the Slavnov- Taylor identity 

m 

k'^'EnwiP^) - mw'^HGip'^) = 0, atp^ = m]j+. (2.8) 

From the above two equations we can get 

RetHG\p2=^2 =0. (2.9) 

Then we obtain the related counterterms of the self-energies as 

5T,HW = R'eS//iy|p2=^2^^ = Re(E{^,^ -h S|^^)|p2=„2^^, 

6T.HG = ReS//G|p2=m2^^ = Re(s{fg S|^C')lp2=m2^^, (2.10) 

where TJ^f^^rHG ™^ ^HVFHG denote the unrenormalized self-energy parts from the fermion 
loop and boson loop diagrams, respectively. The unrenormalized self-energies Tihw and 
T,HG can be calculated from the diagrams in Fig. 2. Due to the on-mass-shell renormal- 
ization condition, the quartic self-energy diagrams in Fig. 2(5) does not contribute to the 
— self-energy. 

Since the incoming photons are unpolarized, it is required to average over their po- 
larization. The sum of the photon polarization vector e^(Ai,pi) and e'^(A2,P2) bas to be 
chosen in such a way that only the physical(transverse) polarization states remain and the 
unphysical one(longitudinal) does not contribute to matrix element] 24]: 



E 6^*(A,p)6^(A,p) = -5^^ + 2^^i^i±^, ip = p^,P2). (2.11) 

A=l,2 ^ 

The renormalized amplitude of the subprocess is expressed as 

= e^(Pi)e''(P2)e^(-^i) {fi9tiuP2X + f29fj,iyPix + f39fi\hu + fiQuxkip, + hK^c^lv\ 

+ hPl(^,iuXa + f7P2^,iuXa + f8h,ikiuP2X + fdhf^hiyPlX + flokfPlP2X(^a,iuf3 

+ fllkfp2P2X^a,iuf3 + fl2kiPlxP2ea,iuf3 + /l3^1 ^li^Pl ea/.A/3 + /uA^^pf ^li.ea/iA/3 

+ hbkiP2klf,eauXp + fwPlP2P2Xefiual3 + fnPlP2PlXei^ual3 + fl8kluPlP2(^fiXal3 

+ fl9h^LPfP2(^uXaf3 + f20kfPiP2eat^f:i-ygiyX + f2lKPlP2<^auf3'ygiiX 

22 

+ f22Kplpleo,x^^9f..) = €''{pi)e''{p2)eH-h) E Z^^- 

i=l 
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where fi{i = 1,22) are the form factors and Ci are their respective Lorentz tensors. The 
u-channel matrix parts ■M^^'^l^, and -M^l^)^ are obtained by the fohowing replacements: 



At (J) {t ^ u,ii ^ u,pi ^ P2), 

bub 
it) 



■M^fl^it ^ u,ii^ U,pi^ P2) 



The amphtude parts of -^(t) (u) ^^"^ -^(tjV) "^^P^^^^^^ contributions from 

the one-loop diagrams Fig.l(a.l-a.l9), the counter term diagrams Fig.l(b.l-b.2) and the 
self-energy bubble diagrams Fig.l(b.3-b.5), respectively. We decompose the amplitude 
parts of M^^-^ and A^^'^"'^ as: 



62 ' 



'64 



'65 



(2.12) 



where 



M 



M 



it) 
61 

62 



M 



(t) 
63 



M 



64 



M 



65 



e^(R)6"fe)e'(-fci) 



2e5V^ 



yWH 



t — m 



H+ 



t — m 



w 



e^^{p^)e^{p2)e\-k,) 



2ie mw 



{t - mfy){t - mjj+] 



'Egw{'^^2 + Cs + 2C4) + t.HW 



2mwJC-3 + {2jC2 -JC-3 + 2^) 

mw 



e^{pi)e''{p2)e\-h) 



(2A + 2/:2 + /:3 + A), 



The explicit expressions of the related unrenormalized self-energies are listed in Appendix. 

The cross section for the subprocess at one-loop order with the unpolarized photon 
collisions can be obtained from 

ri+ 



1 

IGtts^ 



(2.13) 
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In above equation, we define = {m'jj + — s) ± y + — s)^ — Am'jjmyy /2. 
The bar over the sum means average over initial spins. 

With the integrated photon luminosity in the e+e^ collision, the total cross section of 
the process e~^e~ 77 — >^ W~^H~ can be written as: 



cr[s 



dz- 



-77 



o-(77 W^H at s = z s) 



(2.14) 



lEo/y/l dz 

with Eo = mw + ^h+-: ^"^^ \/s(V^) being the e"'"e~(77) center-of-mass energy. dL^^-^jdz 



is defined as: 



—P- = 2z 

dz 



F^/e{x)F^/e{z /x). 

2 I -l^max 



(2.15) 



For the initial unpolarized electrons and laser photon beams, the energy spectrum of 
the back-scattered is given by|25|[26|[27| 



where 



1-X + 



4x 



+ 



4x^ 



l-x C(l-2;) ^2(1 -x) 



(2.16) 



4 8,, ^ 1 

(l----^)ln(l+0 + 7T + 7 



1 



2 e 2(i + e)2 



(2.17) 



me and Eq are the mass and energy of the electron respectively, and ujq is the laser- 
photon energy, x represents the fraction of the energy of the incident electron carried by 
the back-scattered photon. In out evaluation, we choose ujq such that it maximizes the 
backscattered photon energy without spoiling the luminosity through e^e^ pair creation. 
Then we have ^ = 2(1 \/2) ~ 4.8, x^ax ^ 0.83 and D{0 ^ 1.8. 

3 Numerical results and discussions 



In this section, we present some numerical results of the total cross section for the processes 
e~^e~ — > 77 — > W^H^. The SM input parameters are chosen as: mt = 174.3 GeV, 
mz = 91.1882 GeV, mb = 4.4 GeV, sin^ 9w = 0.23117, and oew = 1/128|2|]. 

We take the MSSM parameters in the frame of the mSUGRA scenario. Of the five 
input parameters (mo, mi/2, ^0, tan/? and sign of fi) in this theory, we take mi/2=120 
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GeV, ^0=300 GeV, /x > and set tan (3 to typical values. rriQ is tuned to obtain desired 
mjj± which fall in the experimentally constrained range. 

In the numerical evaluation for the THDM case, we adopted the following MSSM Higgs 
mass relations considering the radiative corrections to fix the neutral Higgs boson masses 
and the mixing angle a p9|. 



tan 2a 



tan 2/?- 



m , 



m| + e/cos(2/?) 



where m\ = m'jj± 



m^z i V ^AZ ~ 4:m^m\cos'^2P — Ae{m\ sin (3 + m\ cos^ 0) 



with 



2)0 F rnj 
V27r^ sin2 /? 



log 



2 



(3.1) 
(3.2) 

(3.3) 



where the left-handed and right-handed scale top-quark masses are determined by formulas 



in a grand unified (GUT) framework |30| 

•f 2 



m,- 



uiq + cos 2/?m^ + TTij 



1^1 + -s^/ cos 2/3m| -I- m^, 



(3.4) 



where at the Plank mass scale we assumed M? = M? = M? = m?,- That is to say in the 
THDM case we take the charged Higgs mass mjj± , niQ and tan (3 as input parameters and 
put their quantities being the same as the corresponding ones in the MSSM case when 
they are allowed in the mSUGRA parameter space, so as to make comparison of the results 
from both models. We noticed also that there are some charged Higgs mass regions in our 
plots which are not covered by the allowed parameter space in the mSUGRA scenario(e.g. 
the dashed lines for the THDM in Fig.3 where tan f3 = 2, 100 GeV < mH± < 300 GeV, 
and mjj± = 200 GeV, 2 < tan /? < 4 in Fig. 4a, et cetera), there we take mo = 220 GeV. 



Actually our evaluation shows that by adopting mSUGRA parameters and Eqs.(p.l|^ |3.4[) 
for the MSSM and THDM models respectively, we can get the corresponding neutral 
Higgs masses and mixing angle a with approximately same quantities in both models. To 
illustrate this point, we list parameter quantities of neutral Higgs masses and mixing angle 
a with tan/3 = 6 used in our calculations for both THDM and MSSM models in table 1. 
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Table 1 The comparison of the corresponding parameter values in the THDM and 

MSSM with tan /3 = 6 



mH± [GeV] 


m/jO [Gey] 


nifjo [GeV\ 


m^o [Gey] 


a 


210 
300 
450 
550 


98.50/96.54 
100.10/97.45 
102.09/97.80 
103.25/97.88 


197.40/197.12 
291.09/290.83 
444.06/443.86 
545.15/544.98 


194.60/193.99 
289.30/289.02 
442.75/442.76 
543.95/544.09 


-0.2547/-0.2644 
-0.2007/-0.2039 
-0.1799/-0.1807 
-0.1750/-0.1753 



Fig. 3 displays the integrated total cross section of W^H^ production at LC with 
y/s = 1 TeV versus the mass of with tan (3 = 2,6, 32, respectively. There we can 
see the sophisticated structures on the MSSM curves, which are induced mainly by the 
threshold effects from Yukawa couplings of charged Higgs boson to quarks and scalar 
quarks, at the vicinities where m}j+ = mi, + (i=l,2) and mu+ = mt + mf, from loop 
diagrams Ip^ . It is obvious that the supersymmetric particles in loop can enhance the 
cross sections of the subprocess substantially by one order. When the ratio of the vacuum 
expectation values has a small value and the charged Higgs boson mass is in the mediate 
range, the cross sections of the subprocess can reach 10^^ femto barn. Since in the case of 
the MSSM we use the mSUGRA senario, while in the THDM we use another parameter 
setting methode mentioned above, some parts of parameter space with small mjj+ and 
tan (3 in the THDM are not allowed in the mSUGRA parameter sapce. That is why we 
see in Fig. 3 that the curves for the MSSM are not plotted at the lower end. 

We present the cross sections of W^H^ productions at LC versus tan (3 in Fig. 4a and 
Fig4.b for a non-SUSY THDM and the MSSM, respectively. Fig.4a is for ^ = 1 TeV, 
the tan (3 is in the range of 2 — 35 and the mass of the charged Higgs boson is taken as 
200 GeV , 400 GeV and 600 GeV , respectively. The Fig. 4a shows that the cross sections 
are reduced with the increment of the charged Hoggs boson mass in both models. In 
Fig. 4a, we can see that when the charged Higgs boson has the mass value about 200 GeV 
and y/s = 1 TeV , the cross sections at LC are in the order of 10~^ femto barn in the MSSM 
THDM and in the range of lO"^ - 10"^ femto barn in the non-SUSY THDM. Fig.4b is 
for \/s = 500 GeV and the mass of the charged Higgs boson is chosen as 200 GeV . The 
two curves correspond to the results for the non-SUSY THDM and the MSSM. All these 
curves in Fig. 4a and Fig. 4b show that the cross section is a obviously effected by tan (3. 
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That is because the couphngs of Higgs bosons to quark(squark) pairs are related to the 
ratio of the vacuum expectation values. We can see from these figures that the dependence 
of the cross section on tan (3 in the MSSM is weaker than in the non-SUSY THDM and 
the maximal enhancement occurs at the low tan/3 range. We compared our results in 
the non-SUSY THDM with those in Ref.|18] and find that the production rates via 77 
collision mode are of the same order as those via e^e~ collision mode at LC. 



4 Conclusion 

In this paper, we study the loop induced VF^'^-associated production of charged Higgs 
bosons via photon-photon fusion at linear colliders in the non-SUSY THDM as well as in 
the MSSM. Numerical analysis of their production rates is carried out at one-loop order 
with some typical parameter sets. With our input parameters, the contributions from 
SUSY sector can lead to an order of magnitude enhancement of the cross section for the 
case of the THDM in some parameter space. Our results demonstrate that this production 
rate via photon-photon collisions at LC can reach about 0.1 femto barn in frame of the 
MSSM at LC when mH± ~ 200 GeV. We find that the H^W^ -associated production 
rate via 77 collisions at LC is quantitatively comparable to that via the electron-positron 
collision mode[p^]. 
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Appendix 

The Feynman rules we used for the couplings can be found in [Q] . We denote the couplings 
concerned in this work as: 
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B+tibj 



V^- ' Pt + V^- ' Pi 



Btib 



BUb, 

rR 



R, 



= V^^o^+Pl + V^^o^+Pr {B+ = H+,G+). 



In our calculation there are several hundred loop diagrams contributing to the sub- 
process 77 ^ W^H^. We find their explicit expressions of one-loop diagrams are too 
complicated and over Icnthy to list here. Therefore, we list only the unrenormalizcd self- 
energies expressions of mixing — W~ and — G~ explicitly in appendix. Their 
fermion and boson parts are expressed as below: 

'iigd 



+ 



32V27r2 

j=l-4 



mtVmbBolp, m, rrib] + {mtVHtb + ■mbVHtb)Bi\p, mt, mb] 
+ Bi\p,m q,m 



HXjXj^ ^3 WXjXj Xi W^Xi' 



327r2 
327r2 



327r2 

327r2 



{Bo [p, m^o , mw] + 2-Bi [p, m/jO , mw] ) sin {a - P) 
■ {Bq \p, tuho , mw] + 2i?i [p, niHO , mw] ) cos (a - /3) 
{Bq [p, mho , mH+] + 2Bi [p, m^o , mH+]) cos (a - /?) 



327r2 

2 



(So [p, m^o , mw] - Bi\p, m^o , mw] ) cos (a - /3) 

{Bo\p, mno , m^^+ ] -I- 2Bi [p, m^^o , ] ) sin (a - /3) 



mwg 
647r2 

E 

ij = l-2 



■(5o[p,"iHO,mM/] - -Bi[p,mj^o,mH/])sin2(a - /3) 

3 3 
^Bo[p,m^.,mg.]F^j-.g.lVg.j. + ^Bi\p,mi^,miJVHi^iy^ 



■J ■ bi' Htj bi Wbitj 

2,d[-Ao[mb]{y}i^,Vl^;, + V§^,Vi;,) +p'B^\p,mt,mb]iV^,b- 

yhib + ymbyctb) + "^t^ob, "^^i (miT/ii^Vc^ti, + ^b ■ 



j=l-4 



—1,2 
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-5o[p,m o,m +]m 

•^7 A-i ^7 



^-HciP^) = {VH+G-hoVH+H-hQBo\p,mho,mH+])/{l6Tr^) 

+ {yG+G-hoVH+G-hoBoip, mhO,mw])/ (IOtt^) 

+ {yG+G-HoVn+G-HoBolp, m^o, mvy])/(167r2) 

- [g^ rn]j+ {Bq [p, niho , mw] - 2Si [p, m^^o , ttivf] + -B21 [p, "ifto , mvt^] ) 



+ 



—dB22[p, rni^o,mw]] sin2(a — /?)} /(1287r^) 

[9^ mjj+{Bo\p,mHO,mw] - 2Bi\p,mH0,mw] + B2i\p, 111^0, mw]) 
—dB22[p, 'mjjo,mw]] sin2(a — /?)} /(1287r^) 
- 3 Bo[p,mi^,mj^jV*-^ V^,^~J{167r^). 

In case of the non-SUSY THDM, we take away the contributions from supersymmetric 
particles, i.e. scalar quarks and charginos/neutralinos, and used dimensional regularization 
(d = 4 — e) instead of dimensional reduction (d=4). As a check of our results, we find from 
our calculation that the amplitude of all the diagrams in Fig.l is UV- finite by themselves. 
The above self-energy expressions for the non-SUSY THDM do not contradict with those 



in Ref.[|l8|. 

In the above expressions we adopt the definitions of one-loop integral functions in 



reference [31|. The numerical calculation of the vector and tensor loop integral functions 
can be traced back to four scalar loop integrals Aq, Bq, Cq, Dq as shown in [|32[. 
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Figure Captions 

Fig.l The Feynman diagrams of the subprocess 77 — > W~^H~ , where the notations 
and 5+ denote neutral particles h^{H^) and charged particles H~^{W~^,G~^), respectively. 

Fig.2 The diagrams contributing to self-energies of mixing i7+ — W~ (Fig.2(l-2)) and 
H+ - G- (Fig.2(3-5)). 

Fig.3 Total cross sections of e'^e" — 77 — W^H^ in mSUGRA scenario as function 
of mH± with = 1 TeV in non-SUSY THDM and the MSSM. 

Fig.4a Total cross sections of e+e~ — 77 — XV^H^ in mSUGRA scenario as function 
of tan/3 with = 1000 GeV in non-SUSY THDM and the MSSM. 

Fig.4b Total cross sections of e"'"e~ — 77 — W^H^ in mSUGRA scenario as function 
of tan^ with ^ = 500 GeV in non-SUSY THDM and the MSSM. 
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